Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 7 December 2009 (MN style file v2.2) 



The Contribution of High Redshift Galaxies to Cosmic 

Reionization: New Results from Deep WFC3 Imaging of 

. . . the Hubble Ultra Deep Field 
a^ ■ 
o ; 

CN ! Andrew J. Bunker^, Stephen Wilkins\ Richard S. EUis^, 
^ '. Daniel Stark ^, Silvio Lorenzoni^, Kuenley Chiu^, Mark Lacy^ 
Q . Matt J. Jarvis^ & Samantha Hickey^ 

^Department of Physics, Denys Wilkinson Building, Keble Road, Oxford 0X1 3RH, U.K. 



O 



email: a.bunkerlOphysics . ox . ac .uk 

2 California Institute of Technology, Mail Stop 169-327, Pasadena, CA 91109, U.S.A. 
^ Institute of Astrophysics, University of Cambridge, Madingley Road, Cambridge, CB3 0HA, U.K. 
NRAO, 520 Edgemont Road, Charlottesville, VA 22903, U.S.A. 



^ j , ^ Centre for Astrophysics, Science & Technology Research Institute, University of Hertfordshire, Hatfield, Herts ALIO 9AB, U.K. 

I Submitted to MNRAS 
O 

C/5 , ABSTRACT 

, We have searched for star-forming galaxies at « 7 — 10 by applying the Lyman-break 

I technique to newly-released J- & //-band images (1.1, 1.25, 1.6 /Ltm) from WFC3 

. on the Hubble Space Telescope. By comparing these images of the Hubble Ultra Deep 

^ ' Field with the ACS z'-band (0.85 /im) images, we identify objects with red colours, 

! {z' — y)AB > 1-3, consistent with the Lyman-a forest absorption at z « 6.7 — 8.8. 

lO , We identify 12 of these z'-drops down to a limiting magnitude Yab < 28.5 (equivalent 

' to a star formation rate of 1.3 MQyr~^ at z = 7.1), all of which are undetected in 

the other ACS filters. We use the WFC3 J-band image to eliminate contaminant low 
Q>>^ ■ mass Galactic stars, which typically have redder colours than z w 7 galaxies. One of 

' our z'-drops is a probably a T-dwarf star. The z « 7 z'-drops appear to have much 

On . bluer spectral slopes than Lyman-break galaxies at lower redshift. Our brightest z'- 

■ drop is not present in the NICMOS J-band image of the same field taken 5 years 

^ I before, and is a possible transient object. From the 10 remaining z « 7 candidates we 

• i-H . determine a lower limit on the star formation rate density of 0.0017 M0yr~^ Mpc~^ 

' for a Salpeter initial mass function, which rises to 0.0025 — 0.004 Mq yr~^ Mpc~'^ after 

^ , correction for luminosity bias. The star formation rate density is a factor of ~ 10 

■ - - ' less than that of Lyman-break galaxies at z = 3 — 4, and is about half the value at 

z « 6. We also present the discovery of 7 F-drop objects with {Y — J)ab > 1-0 and 
Jab < 28.5 which are candidate star-forming galaxies at higher redshifts (z « 8 — 9). 
We find no robust J-drop candidates at z « 10. While based on a single deep field, 
our results suggest that this star formation rate density would produce insufficient 
Lyman continuum photons to reionize the Universe unless the escape fraction of these 
photons is extremely high (/esc > 0.5), and the clumping factor of the Universe is low. 
Even then, we need to invoke a large contribution from galaxies below our detection 
limit (a steep faint end slope). The apparent shortfall in ionizing photons might be 
alleviated if stellar populations at high redshift are low metallicity or have a top-heavy 
initial mass function. 

Key words: galaxies: evolution - galaxies: formation - galaxies: starburst - galaxies: 
high redshift - ultraviolet: galaxies 



1 INTRODUCTION 

In the past decade, the quest to observe the most distant 
galaxies in the Universe has rapidly expanded to the point 



where the discovery of z ~ 6 star-forming galaxies has 
now become routine. Deep imaging surveys with the Hubble 
Space Telescope (HST) and large ground based telescopes 
have revealed hundreds of galaxies at z ~ 6 (Bunker et al. 
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2004, Yan & Windhorst 2004, Bouwcns ct al. 2006, 2007, 
Oesch et al. 2007, Yoshida et al. 2006, McLure et al. 2009). 
These searches typically rely on the Lyman break galaxy 
(LBG) technique pioneered by Steidel and collaborators to 
identify star-forming galaxies at z ~ 3 — 4 (Stcidcl ot al. 
1996, 1999). Narrow-band searches have also proved success- 
ful in isolating the Lyman-a emission line in redshift slices 
between 2 = 3 and z = 7 (e.g., Rhoads et al. 2003, Ota et al. 
2008, Smith & Jarvis 2007, Ouchi et al. 2008), although such 
searches at z > 7 have yet to yield plausible candidates (e.g. 
Willis et al. 2008). Rapid follow-up of Gamma Ray Bursts 
has also detected objects at 2 ~ 6, and most recently one at 
2 = 8.2 (Tanvir et al. 2009, Salvaterra et al. 2009). 

Parallel to these developments in identifying high red- 
shift objects has been the discovery of the onset of the 
Gunn-Peterson (1965) effect. This is the near-total absorp- 
tion of the continuum flux shortward of Lyman-a in sources 
at 2: > 6.3 due to the intcrgalactic medium (IGM) having 
a much larger neutral fraction at high redshift. The Gunn- 
Peterson trough was first discovered in the spectra of SDSS 
quasars (Becker et al. 2001, Fan et al. 2001, 2006). This 
defines the end of the reionization epoch, when the Uni- 
verse transitioned from a neutral IGM. Latest results from 
WMAP indicate the mid-point of reionization may have oc- 
curred at 2 w 11 (Dunkley et al. 2009). The source of nec- 
essary ionizing photons remains an open question: the num- 
ber density of high redshift quasars is insufficient at 2 > 6 
to achieve this (Fan ct al. 2001, Dijkstra et al. 2004). Star- 
forming galcixies at high redshift are another potential driver 
of reionization, but we must first determine their rest-frame 
UV luminosity density to assess whether they are plausi- 
ble sources; the escape fraction of ionizing photons from 
these galaxies, along with the slope of their UV spectra, are 
other important and poorly-constrained factors in determin- 
ing whether star formation is responsible for the ionization 
of the IGM at high redshift. 

Early results on the star formation rate density at z ~ 6 
were conflicting, with some groups claiming little to no evo- 
lution to 2 ~ 3 (Giavalisco et al. 2004, Bouwens et al. 2003) 
while other work suggested that the star formation rate den- 
sity at z ~ 6 was significantly lower than in the well-studied 
LBGs at 2 = 3 — 4 (Stanway et al. 2003). The consensus 
which has now emerged from later studies is that the abun- 
dance of luminous galaxies is substantially less at 2 w 6 than 
at 2 « 3 (Stanway et al. 2003, Bunker et al. 2004, Bouwens 
et al. 2006, Yoshida et al. 2006, McLurc ct al. 2009). If this 
trend continues to fainter systems and higher redshifts, then 
it may prove challenging for star-forming galaxies to provide 
the UV flux needed to fulfill reionisation of the intergalactic 
medium at (e.g. Bunker et al. 2004). Importantly, analysis of 
the faint-end of the luminosity function at high redshift has 
revealed that feeble galaxies contribute an increased fraction 
of the total UV luminosity (Bouwens et al. 2006, 2007, Oesch 
et al. 2008, McLurc et al. 2008), suggesting that the bulk of 
star formation (and hence reionizing photons) likely come 
from lower luminosity galaxies not yet adequately probed 
even in deep surveys. 

Extending this work to the 2 « 7 universe has been 
stunted by small survey areas (from space) and low sensitiv- 
ity (from the ground), limiting current 2 « 7 Lyman-break 
samples to ~ 10 objects (e.g. Bouwens et al. 2008, Oesch 
et al. 2009, Ouchi et al. 2009), none of which have been 



spectroscopically confirmed^. There is evidence of old stel- 
lar populations in 2 ~ 4 — 6 galaxies from measurements 
of the Balmer break in Spitzer/IRAC imaging (Eyles et 
al. 2005,2007; Stark et al. 2007, 2009), which must have 
formed at higher redshift. However, the age-dating and 
mciss-determination of these stellar populations has many 
uncertainties, so searching directly for star formation at red- 
shifts 2 > 7 is critical to measure the evolution of the star 
formation rate density, and address the role of galaxies in 
reionizing the universe. 

The large field of view and enhanced sensitivity of the 
recently-installed Wide Field Camera 3 (WFC3) on HST 
has the potential to make great progress in identifying larger 
samples of 2 > 6 Lyman break galaxies, as it covers an area 
6.5 times that of the previous-generation NICMOS NIC3 
camera, and has better spatial sampling, better sensitivity 
and a filter set better tuned to identifying high-redshift can- 
didates through their colours. In this paper, we present an 
analysis of the recently-obtained WFC3 near-infrared im- 
ages of the Hubble Ultra Deep Field (UDF). We have previ- 
ously used the i'-band and 2'-band ACS images to identify 
LBGs at 2 « 6 through their large i' — z' colours (the i'- 
drops. Bunker et al. 2004). In this paper we use this ACS 
2'-band image in conjunction with the new WFC3 K-band 
to search for galaxies at z « 7, with a spectral break be- 
tween these two filters (the z'-drops). We also anaylse the 
new J and J7-band WFC3 images to eliminate low-redshift 
contaminants of the z'-drop selection through their near- 
infrared colours, and also to look for Lyman-break galaxies 
at higher redshifts (the Y- and J-drops at 2 w 8 & 2 w 10. 

The structure of the paper is as follows. In Section 2 we 
describe the imaging data, the construction of our catalogues 
and our z'-drop selection. In Section 3 we discuss the lumi- 
nosity function of star-forming sources, likely contamination 
by lower-redshift interlopers, and estimate the evolution of 
the UV luminosity function at 2 ~ 7. Our conclusions are 
presented in Section 4. Throughout we adopt the standard 
"concordance" cosmology of Qm = 0.3, Qa = 0.7, and use 
hro = Ho/70kms-^ Mpc'^. All magnitudes are on the AB 
system (Oke & Gunn 1983). 



2 HST/WFC3 OBSERVATIONS 

The near-infrared images of the Hubble Ultra Deep Field 
(UDF) were taken with WFC3 on HST over the period 26 
August - 6 September 2009. The infrared channel of WFC3 
was used, which is a Teledyne 1014 x 1014 pixel HgCdTe de- 
tector, with a field of view of 123 x 136 (a 10-pixcl strip on the 
edge is not illuminated by sky and used for pedestal estima- 
tion). The data form part of the //ST Treasury programme 
GO-11563 (P.I. G. lUingworth). This 192-orbit programme 
will eventually obtain F105W (Y), F125W (J) and F160W 
(H) imaging of the UDF and two deep parallel fields, with 
the UDF observations comprising half of the total orbits. 
The UDF field lies within the Chandra Deep Field South 



^ One galaxy with a spectroscopic rcdsliift oi z = 6.96 was first 
identified through narrow-band imaging for Lyman-a (lye et al. 
2006), but has subsequently also been identified in a survey for 
LBGs, Ouchi et al. 2009). 



Candidate z >7 galaxies from U.ST/WFC3 observations of the UDF 3 



(CDF-S) with the centre of the WFC3 image located at co- 
ordinates RA=03''32™3S4, Decl.=-27°47'00 (J2000), and 
this field is within GOODS-South area (Giavalisco et al. 
2004), surveyed using ACS with the same filters as the ACS 
UDF {b,v,i',z'). The WFC3 field observed falls entirely 
within the 11 arcmin^ of the ACS UDF images. In this pa- 
per, we focus primarily on the F-band WFC3 imaging of the 
UDF obtained over 26-29 August 2009 U.T. This spanned 
18 orbits, split into 2-orbit 'visits' with 2 pointings within 
each orbit. There were small offsets (~< 1 arcscc) between 
the 4 exposures taken in each visit, to allow for sub-pixel 
stepping and to prevent bad pixels repeating. Between vis- 
its there were larger offsets of up to 10 arcscc. The data were 
taken in "MULTIACCUM" mode using SPARSAMPLEIOO, 
which non-destructively reads the array every 100 seconds. 
These repeated non-destructive reads of the infrared array 
allow gradient-fitting to obtain the count rate ("sampling 
up the ramp") and the fiagging and rejection of cosmic ray 
strikes. Each MULTIACCUM exposure comprised 16 reads 
for a total duration of 1403 sec per exposure. 

To guard against the possible contamination of the 
high-redshift «'-drop sample by lower-redshift sources, we 
also analyse the new WFC3 J-band images, as lowor-rcdsliift 
galaxies and Galactic stars will have different Y — J colours 
to high-redshift galaxies even if their z' — Y colours are 
similar. The J-band imaging was taken over 28-30 August 
2009 U.T. over 16 orbits, with each orbit again comprising 
two 1403 sec SPARSAMPLEIOO observations and a similar 
dithering strategy. We also measured the (J — H) colours of 
the z'-drops to determine the UV spectral slopes; the WFC3 
i?-band imaging was taken from 01-05 September 2009 U.T. 
over 28 orbits and 56 SPARSAMPLEIOO exposures. 

2.1 Data Reduction 

The data became public on 9 September 2009, and both raw 
data and the output of the data reduction pipeline for indi- 
vidual frames made available, but not the reduced final com- 
bined images. This pipeline (stsdas . calwf c3 in IRAF) bias- 
and dark-subtracts the images, converts to units of electrons 
by nmltiplication by the gain, and fits the count rate to the 
non-destructive reads of the MULTIACCUM sequence (fiag- 
ging and rejecting cosmic ray strikes in the process). The 
single output frame from the MULTIACCUM sequence is 
fiatfielded, and we chose to work directly with these pipeline 
outputs. Of the 9 visits in the F-band observations, 2 were 
compromised by image persistence - perhaps due to observ- 
ing a bright target shortly before the UDF. These 8 frames 
were excluded from our combined imago in case they intro- 
duced spurious sources. The remaining 28 frames, each of 
1403 sec duration, had residual background removed (which 
might be due to scattered light, variable dark current or 
small errors in the flat-fielding). This was accomplished us- 
ing the XDIMSUM package (M. Dickinson et al.) within IRAF, 
which makes an average background from the disregistered 
(dithered) images, with bright objects masked out so as not 
to bias the background measurement. This average back- 
ground was then subtracted from all the frames, and the 
shifts between the images determined from the centroids of 
several compact objects in the field (the header World Coor- 
dinate System was found to be insufficiently accurate to sat- 
isfactorily register the images). All the frames had the same 



position angle (ORIENT of 128.8 dcg) so this shift-and-add 
was sufficient with no rotation required. The WFC3 point 
spread function at short wavelengths is undersampled by the 
0'.'13 pixels, so we magnified each pixel onto a 3 x 3 pixel grid 
so that we could shift the frames to achieve aligmncnt at the 
sub-pixel level. We combined the background-subtracted, 
rebinned frames with imcombine in IRAF, ignoring pixels 
flagged in a bad pixel mask and rejecting 4cr outliers. We 
survey 4.18arcmin^ in all exposures (39300 sec), and a fur- 
ther 0.67arcmin^ is surveyed with about half of that ex- 
posure time. To make comparison with predicted number 
counts more straightforward, we restrict our analysis to the 
deepest 4.18arcmin^ where the noise is approximately uni- 
form. The J- and f/-band data were reduced in an identical 
manner, except no frames were excluded from the combi- 
nation, as the image persistence effects experienced in 22% 
of the F-band data were not observed in J- and i?-bands. 
Hence the total J-band was 32 exposures (44900 sec), and 
in H it was 56 exposures (78600 sec). 

Usually undersampled imaging data from i/STwould be 
combined through geometric distortion correction and sub- 
pixel registration with the MULTIDRIZZLE software (Koeke- 
mocr et al. 2002). However, at the time of writing, the nec- 
essary distortion files were unavailable for the new WFC3 
instrument. Hence, we determined the mapping between x, y 
pixel space in the WFC3 images, and the celestial coordi- 
nates, by determining the centroids of several hundred ob- 
jects in common between the WFC3 bands and the ACS z'- 
band image of Beckwith et al. (2006) , which had been accu- 
rately astrometrically calibrated. We use the IRAF geomap 
and geotran tasks to fit the distortions with a third order 
polynomial surface, producing residuals of less than 0'.'05. In 
correcting the geometric distortions in the WFC3 images, 
we mapped all the WFC3 Y, J, H images and ACS 6, i' ,z' 
images to the same area, and an output pixel scale of 0'.'06, 
equivalent to 2 x 2 binning of the original drizzled ACS pixels 
(O'.'OS). 

The final frames had units of electrons/sec, and we 
take the standard ACS zeropoints for the UDF images. For 
the WFC3 F105W Y-band, synphot reports an AB magni- 
tude zeropoint of 26.16 (such that a source of this bright- 
ness would have a count rate of 1 electron per second). 
The colour:magnitude diagram (Figure 4) shows reasonable 
agreement between the z'-band and F-band photometry, 
with perhaps a median colour of 0.1 — 0.2 mag (which might 
be due to the average galaxy colour being red, although po- 
tentially the zeropoint might also be uncertain at the ~10 
per cent level). The F125W J-band zeropoint is 26.10, and 
for the F160W iJ-band it is 25.81^ We have checked the 



^ Note added in proof: at the time of writing (December 2009), 
these were still the current zeropoints in use by the MAST 
pipeUne reduction at the Space Telescope Science Institute, using 
the latest calibration files. A web page associated with the WFC 3 
instrument (http: //www. stsci . edu/hst/wf c3/phot_zp_lbn) 
presents slightly different AB zeropoints of 26.27, 26.25 & 25.96 
for F105W, F125W & F160W, which differ by 0.1-0.15 mag from 
the current caUbration files used in the pipeline. Adopting these 
docs not greatly affect tiic selection of Lyman-break galaxies, 
and alters individual star formation rates by Ki 10 per cent, 
comparable to the magnitude uncertainty and less than the 
uncertainty on the conversion of UV flux to star formation rate. 




Figure 1. Postage stamp images of the ^'-drops in the Hubble Ultra Deep Field. The ACS v, i' , and z' images arc shown alongside the 
WFC3 y-band and J-band image for each of our ^'-drops satisfying {z' — Y)j^g > 1.3, Vas < 28.5 and with no detection in the ii-band. 
Each image is 3 arcsec across, with North up and East to the left. 
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Figure 2. The redshift evolution of the {z' — Y) colour of starbursting galaxies, the extrapolated colours of low-redshift galaxies and 
low-mass L and T dwarf stars. The solid, dotted, dashed, and dot-dashed blue lines show the redshift evolution of the colours of 
starburst galaxies with UV continuum spectral slopes f3 S {—2.5, —2.0, —1.5, —1.0, 0.0} (i.e. fx = A'', ranging from no dust reddening 
to ^j^gQQ^ = 4.4 according to the empirical relation of Meurer, Heckman Sz Calzetti 1999). The solid, dotted and dashed red lines 
show the colour of non-evolving low-redshift galaxies (Coleman, Wu & Weedman 1980) as a function of redshift, as three examples of 
foreground objects in our fields (of which elliptical galaxies at z 2.5 would be a contaminant in the z'-drop selection). More discussion 
of contamination by Balmer/4000A break galaxies is given in Wilkins et al. (2009). The yellow and green stars denote the positions of 
L and T dwarf stars respectively (Chiu et al. 2006, Golimowski et al. 2004, Knapp et al. 2004). 



WFC3 //-band photometry against the NICMOS image of 
the same field (Thompson et al. 2005) taken with a very sim- 
iliar F160W filter, and we see good agreement. The NICMOS 
wide-J filter (FllOW) has a significantly difi'erent transmis- 
sion curve to the WFC3 F125W J-band, so to check the new 
J-band photometry we used instead the ESO VLT ISAAC 
images of GOODS-South (Vanzella et al. in prep.). 

In our final combined K-band image, we measure a 
FWHM of 0'.'15 for point sources in the field. As most 
high-redshift galaxies are likely to be barely resolved (e.g.. 
Bunker et al. 2004, Fergusson et al. 2004) we perform pho- 
tometry using fixed apertures of 0'.'6 diameter, and intro- 
duce an aperture correction to account for the flux falling 
outside of the aperture. This correction was determined to 
be 0.2 mag in the F-band and 0.25 mag in J-band and H- 
bands from photometry with larger apertures on bright but 
unsaturated point sources. We note that the //-band im- 
ages display significant Airy diffraction rings around point 
sources. For the ACS images, the better resolution and 
finer pixel sampling require a smaller aperture correction 
of ~ 0.1 mag. All the magnitudes reported in this paper 



have been corrected to approximate total magnitudes (valid 
for compact sources), and we have also corrected for the 
small amount of foreground Galactic extinction toward the 
CDFS using the COBE/DIRBF, & IRAS/ISSA dust maps 
of Schlegel, Finkbeiner & Davis (1998). The optical red- 
dening is E{B ~ V) — 0.009, equivalent to extinctions of 
Afssolp = 0.012, Afiosw = 0.010, Ai25w = 0.008 & 
AiGow = 0.005. 

The geometric transformation and image re-gridding 
produces an output where the noise is highly correlated. 
Hence measuring the standard deviation in blank areas of 
the final registered image will underestimate the noise, so 
to determine the accurate noise properties we consider an 
individual frame of 1400 sec (after fiatfielding, background- 
subtraction and conversion to electrons per second) where 
we measure a standard deviation of ~ 0.02 es~^ in each 
of the combined Y, J and H images. For the final com- 
bined image, a point source of total magnitude Yab = 28.0 
should be detected at 10 cr in our O'.'6-diameter aperture 
(with the aperture correction applied). We confirmed this 
was the case by measuring the noise in a combined version 
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Figure 3. (2' — Y) vs. (Y — J) colour-colour diagram illustrating the differentiation of high-redshift star forming galaxies, low-mass stars 
and lower redshift galaxies. The black circles show the colours of z-dropouts in our selection and their la errors (or lower limits). The 
solid, dotted, dashed, and dot-dashed blue lines show the evolution of the colours of starburst galaxies with UV continuum spectral slopes 
j3 G {—2.5, —2.0, —1.5, —1.0, 0.0} (i.e. f\ = A'' with redshift. The solid, dotted and dashed red lines show the colours of low-redshift 
galaxies (Coleman, Wu & Weedman 1980) as a function of redshift. The yellow and green stars denote the positions of L and T dwarf 
stars respectively (Chiu et al. 2006, Golimowski et al. 2004, Knapp et al. 2004). Note that the J-band photometry of zD3 is suspect, as 
it falls close to the edge of the frame. 



of the data using integer pixel shifts (to avoid smoothing 
the data through interpolation). The faintest sources in our 
catalog (Yab ~ 28.5) should be detected at 6cr. The J- 
and J/-band limits are nearly the same {Jab = 28.0 and 
Hab = 27.9 at 10 cr; and Jab = 28.5 and Hab = 28.4 
at 6 a for point sources measured in a O'.'6-diameter aper- 
ture and with an aperture correction to total magnitude ap- 
plied). At Yab ~ 28.5, a non-detection in the HST z'-band 
(which has a 2 cr sensitivity of z'ab = 30.0) would imply 
a lower limit on the colour of (z' — Y)ab > 1-5 {2a). We 
assess the recoverability and completeness as a function of 
magnitude by adding in fake sources to our images, and de- 
termine the fraction recovered by SExtractor with 3 pixels 
(O'.'lS) of the input source, and with a flux within 0.5 mag. 
We did this by scaling down a bright but unsaturated point 



source and adding this in at different locations on the image 
- most of the high-redshift Lyman-break galaxies show ex- 
tremely compact morphologies well approximated by point 
sources. About 5% of the 4.18arcmin'^ region considered 
was contaminated by foreground galaxies or stars and z'- 
drop candidates in these noisier regions were discounted, 
and our survey volume was corrected for this effect. After 
removing confused regions, the recoverability at our 6 a limit 
{AB ~ 28.5 mag) was 92 per cent, rising to 98 per cent at our 
10a limit {AB ~ 28.0 mag). Hence the completeness correc- 
tions are negligible given our sample size and flux limit. 
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2.2 Construction of Catalogs 

Candidate selection for all objects in the field was performed 
using version 2.5.0 of the SExtractor photometry package 
(Bertin & Arnouts 1996). For the «'-drops, as we are search- 
ing specifically for objects which are only securely detected 
in the WFC3 Y-band, with minimal flux in the ACS im- 
ages, fixed circular apertures 0'.'6 in diameter wore 'trained' 
in the F-image and the identified apertures used to measure 
the flux at the same spatial location in the z'-hand image by 
running SExtractor in dual-image mode. For each waveband 
we used a weight image derived from the exposure map. This 
procedure was repeated for all other ACS and WFC3 filters. 
For object identification, we adopted a limit of at least 5 
contiguous pixels above a threshold of 2cr per pixel (on the 
data drizzled to a scale of Of.' 06 pixel" ^. This cut enabled us 
to detect all significant sources and a rmnibcr of spurious de- 
tections close to the noise limit, or due to diffraction spikes 
of stars or edge effects. Our catalog contains 5000 sources in 
all, 4000 in the region where more than half the frames over- 
lap; many sources are non-unique, split into multiple compo- 
nents due to software deblendin). As high redshift galaxies in 
the rest-UV are known to be compact (e.g., Ferguson et al. 
2004, Bremer et al. 2004, Bunker et al. 2004), we corrected 
the aperture magnitudes to approximate total magnitudes 
with the aperture correction appropriate for that filter. 

2.3 z ^ 7 Candidate Selection 

In order to select z w 7 galaxies, we use the Lyman break 
technique pioneered at z ~ 3 using ground-based telescopes 
by Steidel and co-workers and using HST by Madau ct al. 
(1996). At z ~ 3 — 4 the technique involves the use of three 
filters: one below the Lyman limit (Xrest = 912 A), one in the 
Lyman forest region and a third lougward of the Lyman-a 
line {Xrest = 1216 A). We have shown that by z ~ 6, wc can 
efficiently use only two filters {i'- and z'-band for the z ~ 6 
work), since the integrated optical depth of the Lyman-a 
forest is ^ 1 rendering the shortest-wavelength filter below 
the Lyman limit redundant (e.g., Bunker et al. 2004) . The 
key issue is to work at a sufficiently-high signal-to-noisc ratio 
that drop-outs can be safely identified through detection in 
a single redder band. Here we extend this work to longer 
wavelength, and higher redshift, by using the F-band filter 
on WFC3 to image above the Lyman-a break, and the ACS 
images to find the drop-out objects. In particular, the sharp 
sides of F-band filter (and the SDSS-type filters on ACS), 
coupled the proximity in wavelength of the Y- and z'-bands, 
assist in the clean selection of objects with a sharp spectral 
discontinuity using this photometric redshift technique. In 
Figure 2 we illustrate how a colour cut of (z' — Y)ab > 1.3 
can be effective in selecting sources with z > 6.6. The main 
caveat is that at («' — Y)ab < 2.0 we may be affected by 
contaminants, principally evolved galaxies at z ~ 2.5 (where 
we pick up the Balmer/4000A break), and low-mass cool 
stars in our own galaxy. These L- and T-dwarfs have colours 
of 0.8 < {z' — Y)ab < 1.8, but we may guard against these 
interlopers in our survey by considering the {Y — J) colours 
which are typically much redder than z ~ 7 galaxies with 
similar (z' — Y)ab colours (Figure 3). We emphasize the 
importance in having additional ultrardeep optical imaging 
from ACS (6, v k, i') in the HUDF; any object detected at 



these shorter wavelengths is strongly excluded from being a 
higli-redshift source. 

Our goal is to select objects with colours redder than 
(z' — Y)ab > 1.3 (and preferably redder than (z' — Y)ab = 
1.5, consistent with an LBG at z > 7.0). However, in order to 
ensure we did not exclude gala^xies with colours marginally 
below this (but consistent with the colour cut within the 
magnitude uncertainties), we initially drew up a selection 
with a more liberal colour cut of {z' — Y)ab > 1-0. Imposing 
our selection criteria of {z' — Y)ab > 1-0 and Yab < 28.5 on 
the catalog resulted in 110 objects, 24 of which were clearly 
spurious (mainly diffraction spikes from stars) . A further 40 
had intermediate colours of 1.0 < (z' — Y)ab < 1-3, and 
only 2 of these were within 1 cr of (z' — Y)ab ~ 1-5 - so this 
sub-sample is probably low redshift. We note that 5 objects 
in our intermidate colour sample, 1.0 < (z' — Y)ab < 1-3, 
are undetected in the b, v, i' ACS images (with 6as < 30.3, 
vab < 30.7 and i'AB < 30.6 at 2 a) and may well be 
i'-drop galaxies at the high-redshift end of this selection 
z ~ 6.3 — 6.5, where the break occurs in the z'-band (leading 
to red i' — z' and z' — Y colours); indeed, two of these sources 
appears in the Yan & Windhorst (2004) catalog of faint i!- 
drops in the UDF (their sources #75 & 85, which also appear 
in the catalog of Bouwens et al. 2006) . The other 3 sources 
with intermediate (z' — F ) colours and non-detections in the 
h,v,i' have coordinates 03:32:41.83 -27:46:11.3, 03:32:37.44 
-27:46:51.3 & 03:32:39.58 -27:46:56.5, and magnitudes 
Yab = 28.32 ± 0.09, 27.75 ± 0.06 & 27.87 ± 0.06. We ex- 
clude these intermediate-colour objects from our z > 6.5 
selection. The remaining 46 objects with {z' — Y)ab > 1-3 
were then analysed further. We demanded that these sources 
be undetected in the ACS optical images, as the v-band falls 
below the Lyman limit at the redshifts of interest. Twenty 
galaxies were detected at > 2 a in v-band and i'-band (with 
vab < 30.7 and i'ab < 30.6) and were eliminated as low- 
redshift sources. A further 3 were detected in i' but not in 
V and arc probably not z ~ 7 galaxies. Of the remaining 23 
objects, detailed examination eliminated 11 as probably spu- 
rious (due to fainter diffraction spikes, noisy areas of the im- 
age, or extended outer regions of clearly low- redshift galaxies 
which had been identified by SExtractor as drop-outs due 
to the different PSFs between WFC3 and ACS). This left 
our core sample of 12 z'-drops. We consider whether any of 
these 12 might still be low redshift objects in Section 3.1. 



3 ANALYSIS 

3.1 Discussion of Candidate z m 7 galaxies 

We have 12 candidates meeting our selection criteria of 
(z' — Y)ab > 1.3 and Yab < 28.5, after removing clearly 
spurious sources and those detected at shorter wavelengths 
in the ACS bands. We now consider whether any of these 
12 z'-drops are plausibly L- or T-dwarf Galactic stars. Since 
the resolution of WFC3 is less good than ACS, and these ob- 
jects are typically undetected by ACS, morphological cuts 
(i.e. eliminating apparently unresolved point sources) are 
not appropriate as they would also remove compact galax- 
ies. Hence we use the z'YJ colours to differentiate potential 
stars from the high-redshift galaxy tracks in colour:colour 
space (Figure 3). Of our 12 candidates, object zDll has 
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ity of our candidates, except those brighter sources previ- 
ously identified by Bouwens et al. (2004,2008), are unde- 
tected in the NICMOS image (which does not go as deep 
as the new WFC3 image). However, our brightest z'-drop 
candidate (zDO), should have easily been detected, but it is 
not present in either the NICMOS J ox H images (which 
reach Jab = 26.7 and Hab = 26.2 at 5cr), despite falling 
within the field of view and having Jab{F125W) = 25.6. 
This compact source is well-detected in all of the Y-, J- and 
_ff-band WFC3 images, so is clearly a real object rather than 
an artifact. It appears possible that this may be a transient 
object, perhaps a supernova. We note that the Yan & Wind- 
horst (2004) i'-band drop-out objects 56a& b arc the closest 
well-detected ^'-barid objects (l'.'2 away), but it is unclear if 
there is a physical connection between these z f» 6 sources 
and the transient observed with WFC3. 



20 22 24 26 28 30 
Y(AB) 

Figure 4. The (2' — Y)ab colours of all objects in our catalog 
(dots), as a function of Y-band magnitude. The vertical dotted 
line indicates a lOir detection (Yab = 28.0). Our primary z'- 
drop selection region where (z' — Y)ab > 1-3 and Yab < 28.5 is 
marked by the dashed lines - candidates lie in the upper left, and 
those z'-drops which are undetected at short wavelengths (absent 
from the ACS i)-band) are marked as crosses, or arrows denoting 
2 <T lower limits on their colours if they are undetected in z'-band. 



colours consistent with being a T-dwarf, although we stress 
that it is undetected in the z'-band and the {z' — Y)ab > 1-5 
limit means it could still lie at high redshift. However, to be 
conservative, wo regard this object as a probable contami- 
nant. We note that, as our faintest z'-drop, if we were instead 
to retain it in our analysis it would have negligible impact 
(raising the integrated counts to Yab = 28.5 by 10 per cent, 
and increasing the integrated star formation rate density 
by 5 per cent). All the other candidates have photometry 
distinct from that expected from low-mass stars, with the 
possible exception of zD9 which lies > 2cr away in colour 
space from the most extreme stellar colours. Both zD9 and 
zDll are compact, but sufficiently faint in the near-infrared 
that star-galaxy separation is inaccurate; zD9 appears to 
have a slight Ecist-West elongation, which would rule out 
the stellar interpretation. 

Three of our candidates appear in previous studies 
(zDl, zD2 and zD3); each of these were identified previously 
as 2;'-band drop outs in Bouwens et al. (2004) (UDF-640- 
1417, UDF-983-964, UDF-387-1125), based on the shallower 
NICMOS imaging in FllOW (J- wide) and F160W (H), 
along with the same UDF ACS z'-band as used here. The 
NICMOS photometry from Bouwens et al. (2004) is broadly 
consistent with our WFC3 measurements in Y F105W, J 
F125W and H, although we note the filter bandpass shape 
of NICMOS FllOW J- wide is very different from the WFC3 
filters. Source zD2 is also presented in Oesch et al. (2009), 
and source zDl is red object #3 in Yan & Windhorst (2004). 
We have also inspected the old NICMOS images of the 
UDF (Thompson et al. 2005), and note that the major- 



3.2 UV luminosity function at z ~ 7 

Wo use the observed surface density of our z'-drops in the 
UDF (using a catalog cleaned of spiuious sources, and cor- 
recting the area surveyed for that masked by foreground 
objects, w 5 per cent) to compare with previous estimates 
of the rest-UV luminosity function at high redshift. Our 
final z --^ 7 candidate list comprised 10 objects, with the 
likely star (zDll) and probable transient (zDO) removed^. 
We used the Schechter luminosity function with different 
values of faint end slope (a), characteristic number density 
{(f>*) and luminosity (L*) to predict the observed number 
of sources as a function of magnitude, which would obey 
our (z' — Y) colour cut. We assume a simple model spec- 
trum for these high-redshift star-forming galaxies, with a 
99 percent flux decrement below Lyman-a (i.e. Da = 0.99) 
due to the Lyman-a forest during the Gunn-Peterson era. 
We adopt a power law at longer wavelengths, with fx oc X^. 
Measurements of the spectral slopes of i'-drops at z ~ 6 indi- 
cate rather blue colours compared to lower redshift Lyman- 
brcak galaxies (Stanway, McMahon & Bunker 2005), with 
P ^ —2 —2.2 (i.e. approximately flat in f^) and indeed 
the z'-drops at z ~ 7 also seem to exhibit very blue spec- 
tral slopes (Figure 6). Such blue colours might indicate lit- 
tle or no dust extinction for a constant star formation rate, 
according to the empiral relation 0,-^^^^^ = 4.43 -|- 1.99/3 
(Meurer, Heckman & Calzetti 1999), but could also be pro- 
duced through a top-heavy IMF or low metallicity. For a 
given limiting magnitude, our code considers redshift slices 
over which Lyman-broak galaxies would satisfy the colour 
cut, and computes the equivalent limiting luminosity (in 
terms of L*) and comoving volume probed. These are then 
integrated over all redshifts to produce predicted number 
counts for comparison with our observations. 

^ Note added in proof: an independent study of the HUDF- 
WFC 3 images has been posted to the preprint server (Oesch 
et al., submitted arXiv:0909.1806), which uses slightly bluer 
colour cuts and goes to lower signal-to-noise. All 10 of the z'- 
drops presented here also appear in the Oesch et al. sample, 
along with another four which are too blue for our selection (at 
0.8 < {z' — Y)ab < 1-0, and presumably lie at slightly lower 
redshift than our sub-sample) and another two which are fainter 
than our Yab < 28.5 cut. 
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We consider three luminosity functions. The well- 
studied Z7-band drop-out Lyman-break galaxies at 2; ~ 3 
have a rest-framo UV luminosity function with M* = = 

1500A 

-21.1 (equivalent to SFR*=3 = ISMgyr"^), 0*=3 = 
0.00138 Mpc"^ and a faint-end slope of a = -1.6 (Steidel 
et al. 1996, 1999). Our a'-drops exhibit spectral slopes close 
to /3 = —2.0 (i.e., flat in /^) meaning that the ^-correction 
to our observed rest-frame wavelength of ~ 1350 A is negli- 
gible. We model the predicted number counts in the instance 
of no evolution in the properties of the star-forming popu- 
lation from 2; ~ 7 to a ~ 3. We also consider the luminosity 
function derived from a survey of 506 i'-drops at 2 ~ 6 by 
Bouwens et al. (2006) , who conclude that the faint end slope 
is steeper (a = —1.73), with strong evolution in the lumi- 
nosity (M*=g = —20.25, a factor of two fainter than M*=3) 
and little evolution in the characteristic number density 
= 0.00202 Mpc"'^). Finally, wo consider the effects of 
the suggested trend of fainter M* at higher redshifts, adopt- 
ing the parameterization of M^jy = —21.0 -|- 0.36 x {z — 3.8) 
suggested by Bouwens et al. (2009) on the basis of the 
paucity of a'-drops in HST imaging with NICMOS (more 
shallow than the WFC3 data analysed here), for a predicted 
value of M*=7 = -19.8 (with = 0.0011 Mpc"^ and 

a = -1.73). 

In Figure 5 we compare the number counts of galaxies 
with colours (2' — Y)ab > 1-3 to the predictions based on 
the model luminosity functions at 2 = 3, 6, 7. As can be seen, 
our observed cumulative surface density is 1.9±0.7 arcmin"^ 
brighter than Yab = 28.0 (our 10 cr limit. This is well below 
that predicted assuming the 2 = 3 UV luminosity function 
describes galaxies at 2 ~ 7 (it would predict 6.9 arcmin"'^)*. 
However, our observed surface density is comparable to 
the luminosity function estimates assuming the 2 = 6 and 
2 = 7 models (3.7arcmin~^ and l.Oarcmin"^ respectively). 
Pushing to fainter magnitudes (where our completeness de- 
creases), our surface density of 2'-drops with Yab < 28.5 
is 2.4 ± 0.7arcmin~^, which is significantly below both the 
2 = 3 and 2 = 6 predictions of 11.5 & 7.3 arcmin"^, and 
consistent with the 2 = 7 model prediction of 2.2arcmin~^. 
The volume we probe is relatively small (~ 10^ Mpc"^ co- 
moving), and given the observed source counts, we estimate 
a cosmic variance of 35% assuming a pencil beam survey 
geometry with area 4.18 arcmin^ and a mean redshift of 7.0 
with A2 = 1.0 (Trent! & Stiavolh 2008). This uncertainty 
due to Cosmic Variance is comparable to the shot noise on 
the small number of galaxies in our sample (l/VTO). 

Hence we strongly rule out the simple scenario of no 
evolution over the range 2 = 7 — 3 as the observed counts 
arc 3 5 times too low, as we also established at 2 6 with 
the i'-drops (Stanway, Bunker & McMahon 2003; Bunker et 
al. 2004). At the bright end, this interpretation agrees with 
the low surface density of z'-band drop-outs found in the 
HAWK-I imaging of GOODS-South presented in Hickey et 



Recently, Rcddy &z Stcidcl (2009) have re-fit a steeper lumi- 
nosity function to the z f» 3 Lyman-break population, with 
a = -1.73, (p* = 0.00171 and M* = -20.97. We have 

recomputed the expected surface density of objects obeying our 
colour cuts for this luminosity function at 2 f» 7, and the change 
is small: the 17% more t Yab < 28 and 25 per cent higher at 
Yab < 28.5. 



CD 




E 



3 

z r , ,/ / I 

25 26 27 28 29 

mag lim (Fab) 

Figure 5. The observed cumulative surface density of (z' — 
Y)ab > 1-3 objects as a function of limiting magnitude in Y- 
band. Overplotted are three rest-UV luminosity functions for 
high-redshift Lyman-breaJj galaxy populations - the solid curve 
is for the f7-band drop-outs at 2 ~ 3 (Steidel et al. 1996, 1999); 
and the dotted curve is for the i'-drops at z ~ 6 (Bouwens et 
al. 2006). The dashed curve represents a prediction at z = 7 in 
a scenario where M* fades linearly with (l-|-z) (Bouwens et al. 
2009). 

al. (2009) of w 0.01 - 0.02 a.rcmin''^ at Yab < 26.0.We also 
find evidence that the 2 ~ 6 luminosity function of Bouwens 

et al. (2006) over-predicts the number of 2'-drops at 2 = 7 
by a factor of 2 — 3, suggesting evolution in the star forming 
population from 2 = 7 — 6. 

3.3 The Star Formation Rate Density at 2 w 7 and 
Implications for Reionization 

We can use the observed F-band magnitudes of objects in 
the 2'-drop sample to estimate their star formation rate from 
the rest-frame UV luminosity density. In the absence of dust 
obscuration, the relation between the flux density in the rest- 
UV around ~ 1500 A and the star formation rate (SFR in 
MQyr-^) is given by Lvv = 8 x W'^'^SFRergs s'^ Hz''^ 
from Madau, Pozzetti & Dickinson (1998) for a Salpeter 
(1955) stellar initial mass function (IMF) with 0.1 M© < 
M* < 125 Mq. This is comparable to the relation derived 
from the models of Leitherer & Heckman (1995) and Ken- 
nicutt (1998). However, if a Scalo (1986) IMF is used, the 
inferred star formation rates will be a factor of « 2.5 higher 
for a similar mass range. In the absence of a spectroscopic 
redshift, wo assume that these lie at the predicted aver- 
age redshift for galaxies obeying our colour cuts and with 
Yab < 28.5. For the luminosity functions considered, the 
predicted mean redshift is around 2 = 7.1 for a spectral 
slope l3 ~ —2.0. At the limit of our catalog (Yab = 28.5), 
the inferred star formation rate for a galaxy at 2 = 7.1 is 
1.3MQyr"\ equivalent to 0.1 L' for M^y = -21.1 (as at 
2 = 3) or 0.2 L* for M^y = -20.24 (for the 2 = 6 luminosity 
function of Bouwens et al. 2007). The total star formation 
rate in our 10 candidate 2 ~ 7 galaxies is 29.6 Mq yr^^ as- 
suming they are at the predicted mean redshift. The colour 
selection of (2' — Y) is effective over the redshift range 6.7 < 
2 < 8.8 (Figure 2), and so the 4.2arcmin^ surveyed in the 
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ID 


RA & Declination J2000 


Yae 










Jab 






Hab 




(z' - Y)ab 




-J)ab 




zDO" 




53 


OV. /I '7. Qi? 


.0 


26.43 


1 

± 


0. 


,03 


29.95 ± 0.62 


25, 


,62 ± 0, 


,01 


25, 


.77 


1 

± 


0, 


,02 


3.52 ± 0.62 


0.81 


± 


0.03 


[8.82] 


zDl'' 


03:32:42, 


,56 


-27:46:56, 


,6 


26.53 


± 


0. 


,03 


28.05 ± 0.11 


26, 


,23 ± 0, 


,02 


26, 


.26 


± 


0, 


,02 


1.53 ± 0.11 


0.29 


± 


0.04 


8.05 


zD2'' 


03:32:38. 


,81 


-27:47:07, 


.2 


27.35 


± 


0. 


,06 


>30.0 (20-) 


26, 


,74 ± 0, 


,03 


26, 


.66 


± 


0, 


,03 


>2.65 (2o-) 


0.60 


± 


0.07 


3.79 


ZDS'" 


03:32:42. 


.57 


-27:47:31, 


.5 


27.46 


± 


0. 


,06 


29.35 ± 0.36 




27.12= 




27 


.19 


± 


0, 


,05 


1.89 ± 0.36 


0.34 


± 


0.07 


3.42 


zD4 


03:32:39. 


.55 


-27:47:17 


.5 


27.84 


± 


0. 


,09 


29.68 ± 0.48 


27 


.34 ± 0. 


,05 


27 


.15 


± 


0, 


.05 


1.84 ± 0.49 


0.50 


± 


0.10 


2.41 


zD5 


03:32:43. 


.14 


-27:46:28, 


.5 


27.95 


± 


0. 


,11 


>30.0 (2a) 


27, 


.74 ± 0. 


,07 


27, 


.65 


± 


0, 


.08 


>2.05 (2(t) 


0.20 


± 


0.13 


2.19 


zD6 


03:32:36. 


.38 


-27:47:16, 


.2 


27.95 




0. 


,11 


29.32 ± 0.35 


27, 


.83 ± 0. 


,08 


28, 


.11 


± 


0, 


.12 


1.37 ± 0.37 


0.12 


± 


0.14 


2.18 


zD7 


03:32:44. 


.70 


-27:46:44, 


.3 


27.95 


± 


0. 


,11 


>30.0 (2(t) 


27, 


.08 ± 0, 


,04 


27, 


.32 


± 


0, 


.06 


>2.05 {2(t) 


0.87 


± 


0.12 


2.18 


zD8 


03:32:40, 


,57 


-27:46:43, 


.6 


27.98 




0. 


11 


29.45 ± 0.39 


27 


,97 ± 0, 


,09 


28, 


.33 


± 


0, 


.15 


1.48 ± 0.41 


0.01 


± 


0.14 


2.13 


zD9= 


03:32:37, 


,21 


-27:48:06, 


.2 


28.17 


± 


0. 


13 


>30.0 (2(7) 


27 


,63 ± 0, 


,07 


27 


.70 


d= 


0, 


.08 


>1.83 (2o-) 


0.53 


± 


0.15 


1.78 


zDlO 


03:32:39, 


,73 


-27:46:21, 


.3 


28.40 


± 


0, 


,16 


>30.0 (2c7) 


28, 


,29 ± 0, 


,12 


28, 


.53 


lb 


0, 


.18 


>1.6 (2(t) 


0.11 


± 


0.20 


1.44 


zDll'' 


03:32:38, 


,36 


-27:46:11, 


.9 


28.45 


± 


0, 


,16 


>30.0 (2cr) 


27, 


,86 ± 0, 


,08 


28, 


.33 


± 


0, 


.15 


>1.55 (2cr) 


0.58 


± 


0.18 


[1.38] 



t The star formation rate (SFR) is in units of Mq yr ^ based on a relation between the rest-UV continuum density of 
• = 8 X lO^^ergs"^ / SFR, where we assume that the objects are at the predicted average redshift of the z'-drops {z = 7.1). 
Where the SFR is in brackets, these objects are unlikely to be at z Ri 7, but the equivalent star formation rate for their Y-band 



magnitudes is listed but not used in our calculation of the stax formation rate density. 

Table 1. z-band drop out candidates meeting our selection criteria with Y < 28.5 AB mag. ° likely a SN. identified as z-band drop 
outs in previous studies. possibly a star based on their Yios — J125 colour. likely to be a star based on their Yios — J125 colour. " 
The J-band photometry of this source (zD3) is compromised as it falls on the edge of the frame. 



ID 


RA & Declination J2000 


Jab 


Yab 


Hab 


(Y - J)ab 


(J - H)ab 


Vhi 1 arcsec 


YDl 


03:32:42.88 


-27:46:34.5 


27.70 ± 0.08 


28.89 ± 0.24 


27.98 ± 0.10 


1.19 ± 0.25 


-0.28 ± 0.13 


0.38 


YD2 


03:32:37.80 


-27:46:00.1 


27.88 ± 0.10 


30.12 ± 0.73 


28.07 ± 0.11 


2.24 ± 0.73 


-0.19 ± 0.14 


0.14 


YD3 


03:32:38.14 


-27:45:54.0 


28.07 ± 0.11 


29.77 ± 0.53 


28.07 ± 0.11 


1.70 ± 0.54 


0.00 ± 0.15 


0.16 


YD4 


03:32:33.13 


-27:46:54.4 


28.11 ± 0.11 


29.74 ± 0.66 


28.9 ± 0.23 


1.63 ± 0.67 


-0.79 ± 0.26 


0.11 


YD5 


03:32:35.85 


-27:47:17.1 


28.38 ± 0.14 


29.44 ± 0.39 


27.98 ± 0.10 


1.06 ± 0.42 


0.41 ± 0.17 


0.10 


YD6 


03:32:40.10 


-27:47:18.8 


28.40 ± 0.16 


29.99 ± 0.67 


28.07 ± 0.19 


1.59 ± 0.69 


-0.31 ± 0.25 


0.11 


YD7 


03:32:37.63 


-27:46:01.5 


28.44 ± 0.16 


29.57 ± 0.43 


28.61 ± 017 


1.13 ± 0.46 


-0.17 ± 0.24 


0.10 



Table 2. Coordinates and magnitudes of the Y-band drop out candidates meeting our selection criteria with (Y — J)ab > 1-0, J < 28.5 
AB mag and no detection at > 2 cr in the ACS images. We also present the measured half-light radius (for an unresolved source, 
Thi = 0'.'08). 



deepest region of the WFC3 Y-imago should colour-select 
a'-drop galaxies within a cornoving volume of 17700 Mpc'^. 
So we calculate the star formation rate density from the 
a'-drops in our sample to be 0.0017 yr"^ Mpc~^. How- 
ever, this should bo regarded as a conservative lower limit on 
the star formation rate density, because the redshift range 
is not surveyed with uniform sensitivity to UV luminosity 
(and hence star formation rate). We have a strong luminosity 
bias towards the lower end of the redshift range, due to the 
effects of increasing luminosity distance with redshift and 
also the Lyman-a break obscuring an increasing fraction of 
the filter bandwidth. Using an "effective volume" approach 
(e.g., Steidel et al. 1999) accounts for this luminosity bias, 
and increases the inferred total star formation rate density 
to 0.0034 Afoyr^^Mpc^'"*, integrating down to 0.2 Luy of 
the L* value at 2; = 3 (i.e. Muv = —19.35) and assuming 
a UV spectral slope of /3 = —2.0. If instead we integrate 
down to 0.1LJ=3 (0.2L*=e) and the total star formation 
rate density is 0.004 M© yr~^ Mpc~^. These star formation 
rate densities arc a factor of ^ 10 lower than at ^ 3 — 4, 
and even a factor of 1.5 — 3 below that at 2; « 6 (Bunker et 
al. 2004; Bouwens et al. 2006). 

We can compare our measured UV luminosity density at 
« « 7 (quoted above as a corresponding star formation rate) 
with that required to ionize the Universe at this redshift. 
Madau, Haardt & Rees (1999) give the density of star format 



tion required for reionization (assuming the same Salpeter 
IMF as used in this paper): 

0.005 Mq yr"^ Mpc"^ fl + z\^ ( Q.bh^mV (C\ 

^'^^'^ 1::^ V^) \0^) I5J 

We have updated equation 27 of Madau, Haardt &: Rees 
(1999) for a more recent concordance cosmology estimate of 
the baryon density from Spergel et al. (2003). 

C is the clumping factor of neutral hydrogen, C = 
{pui} {ptii)~^ ■ Early simulations suggested C « 30 (Gnedin 
& Ostriker 1997), but more recent work including the effects 
of reheating implies a lower concentration factor of C w 5 
(Pawlik et al. 2009). The escape fraction of ionizing pho- 
tons (/esc) for high-rcdshift galaxies is highly uncertain (e.g., 
Steidel, Pettini & Adelbergcr 2001, Shapley et al. 2006), but 
even if we take /esc = 1 (no absorption by H i) and a very 
low clumping factor, this estimate of the star formation den- 
sity required is a factor of 1.5-2 higher than our measured 
star formation density at 2; « 7 from z'-drop galaxies in the 
UDF. For faint end slopes of a — 1.8 —1.6 galaxies with 
L > 0.2 L* account for 24 — 44% of the total luminosity (if 
there is no low-luminosity cut-off for the Schechter function) , 
so even with a steep faint-end slope at 2; « 7 we still fall 
short of the required density of Lyman continuum photons 
required to reionize the Universe, unless the escape fraction 
is implausibly high (/esc > 0.5) and/or the faint end slope is 
a < —1.9 (much steeper than observed at 2: = — 6). How- 
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ever, the assumption of a solar metallicity Salpetcr IMF may 
be flawed: the colours of a ~ 6 i'-band drop-outs arc very 
blue (Stanway, McMahon & Bunker 2005), with /3 < -2, 
and the new WFC3 J- and H-hand images show that the 
z ~ 7 z'-drops also have blue colours on average (Figure 6). 
A slope of /3 < — 2 is bluer that for continuous star forma- 
tion with a Salpeter IMF, even if there is no dust reddening. 
Such blue slopes could be explained through low metallicity, 
or a top-heavy IMF, which can produce between 3 and 10 
times as many ionizing photons for the same 1500 A UV lu- 
minosity (Schaerer 2003 - see also Stiavelli, Fall & Panagia 
2004) . Alternatively, we may be seeing galaxies at the onset 
of star formation, or with a rising star formation rate (Vorma 
et al. 2007), which would also lead us to underestimate the 
true star formation rate from the rest-UV luminosity. We 
explore the implications of the blue UV spectral slopes in 
z > 6 galaxies in a forthcoming paper (Wilkins et al. in 
prep.). 



3.4 Star-Forming Galaxies at z ~ 8 and Beyond 

The availability of deep longer- wavelength WFC3 images in 
addition to the K-band enables us to push the Lyman break 
technique to higher redshifts still, searching for z « 8 objects 
where the Lyman-a decrement falls below the J-band ("Y- 
drops"), and z « 10 sources where the break falls between 
the J- and _ff-filters ("J-drops"). Again using SExtractor, 
but training in the J-images as our detection band, we un- 
cover 7 objects with {Y — J)ab > 1-0, Jab < 28.5 and no 
ACS detection (Figure 7 and Table 2). These F-drops are 
potential z ~ 8 star-forming galaxies, and a discussion of 
the selection and properties of these z ~ 8 candidates and 
the possible contamination fraction will appear in a forth- 
coming paper (Lorenzoni et al. in prep.)^. Using H as our 
detection band, we find no robust (J — Hab) > 1.0 objects 
down to Hab < 28.5 which are undetected by ACS. The 
absence of robust J-drops hints at further evolution of the 
UV luminosity function at z S> 7. 



4 CONCLUSIONS 

We have searched for star-forming galaxies at z « 7 by ap- 
plying the Lyman- break technique to newly-released 1.1 /xm 
F-band images from WFC3 on the Hubble Space Telescope. 
By comparing these images of the Hubble Ultra Deep Field 
with the ACS z'-band (0.85 /xm) images, we select objects 
with red colours, (z' — Y)ab > 1-3, consistent with the 
Lyman-a forest absorption at z w 6.7 — 8.8. We identify 
12 of these z'-drops down to a limiting magnitude Yab < 
28.5 (equivalent to a star formation rate of 1.3 Mq yr^^ at 
z = 7.1), all of which are undetected in the other ACS fil- 
ters, consistent with the Lyman limit break at z w 7. We 
also analyse the new 1.25 /im J-band and 1.6 /xm Jif-band 
WFC3 images to measure the near-infrared colours of our 



'' Note added in proof: an independent analysis of y-drops in the 
WFC 3-HUDF has been posted to the preprint server (Bouwcns et 
al. submitted, arXiv:0909.1803); our sample contains all but one 
of the Y-drops presented in Bouwens et al. , with UDFy-43086276 
being fainter than our Jab < 28.5 cut. As with our analysis, 
Bouwens et al. find no robust J-drop candidate at z ~ 10. 
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Figure 6. The near-infrared colours of the z'-drops in the Hubble 
Ultra Deep Field. Many have very blue (J — H) colours. The 
dotted horizontal lines show the (J — H) colours produced by 
spectral slopes of /3 = —1, —2, —3, and the dashed vertical lines 
show the redshifts which would produce the {z' — Y) colours, 
assuming no emission line contamination from Lyman-a and a 
spectral slope of /3 = —2.0. The triangle is the probably transient 
source (zDO) and the square is the likely T-dwarf star (zDll). 



z'-drops, and eliminate low mass Galactic stars which typ- 
ically have redder colours than z « 7 galaxies. One of our 
z'-drops is a probably T-dwarf star. Our brightest z'-drop 
is not present in the NICMOS J-band image of the same 
field taken 5 years before, despite being well above the flux 
limit of this less sensitive image. This is a possible transient 
object, and we exclude this and the probable Galactic star. 

From the 10 remaining 2 ~ 7 candidates we deter- 
mine a lower limit on the star formation rate density of 
0.0017 M© yr~^ Mpc~^ for a Salpeter initial mass function, 
which rises to 0.0025 - 0.0034 Mq yr"^ Mpc~^ after correc- 
tion for luminosity bias. The star formation rate density is 
a factor of ~ 10 less than that of Lymari-break galaxies at 
z = 3 — 4, and is about half the value at z « 6. This star for- 
mation rate density would produce insufficient Lyman con- 
tinuum photons to reionize the Universe unless the escape 
fraction of these photons is extremely high {fcsc > 0.5) and 
the clumping factor of the Universe is low. Even then, we 
need to invoke a large contribution from galaxies below our 
detection fimit (i.e., a steep faint end slope). The apparent 
shortfall in ionizing photons might be alleviated if stellar 
populations at high redshift are low metallicity or have a 
top-heavy initial mass function. 

We have also inspected the YJH images for the pres- 
ence of Lyman-break galaxies at higher redshifts. We iden- 
tify 7 objects with (Y — J)ab > 1-0 down to Jab < 28.5; 
these Y-drops are candidate z « 8 star-forming galaxies. 
A similar selection based on the i/-band did not find any 
J-drops with (J — H)ab > 10 and Hab < 28.5 - all such 
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Figure 7. Postage stamp images of the Y-drops in the Hubble Ultra Deep Field. The ACS v, i', and z' images are shown alongside the 
WFC3 Y, J and //-bands for each object satisfying (Y — J)ab > l-Oi Jab < 28.5 and with no detection in the ACS images. Each image 
is 3 arcsec across, with North up and East to the left. 



objects were detected in the ACS images (below the Lyman- 
limit at 2 ~ 10), and ruling out a high-redshift interpreta- 
tion. 
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